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Abstract 

A neural net architecture is proposed to estimate, in real-time, the fatigue life of mechanical 
components, as part of the Intelligent Control Systems for Reusable Rocket Engines. Arbitrary 
component loading values were used as input to train a two hidden-layer feedforward neural net 
to estimate component fatigue damage. The ability of the net to learn, based on a local strain 
approach, the mapping between load sequence and fatigue damage has been demonstrated for 
a uniaxial specimen. Because of its demonstrated performance, the neural computation may 
be extended to complex cases where the loads are biaxial or triaxial, and the geometry of 
the component is complex (e.g. turbopump blades). The generality of the approach is such 
that load/damage mappings can be directly extracted from experimental data without requiring 
any knowledge of the stress/strain profile of the component. In addition, the parallel network 
architecture allows real-time life calculations even for high frequency vibrations. Owing to its 
distributed nature, the neural implementation will be robust and reliable, enabling its use in 
hostile environments such as rocket engines. This neural net estimator of fatigue life is seen 
as the enabling technology to achieve component life prognosis, and therefore would be an 
important part of the life extending control for reusable rocket engines. 


1 Introduction. 

Intelligent Control Systems (ICS) for Reusable Rocket Engines (RRE), such as the Space Shuttle 
Main Engine (SSME), require expert system based decision-making within the response times of 
the RRE [1]. For a given level of knowledge of the RRE, and a given amount of computational 
power allocated to the expert system, accurate decision-making requires a complete and timely 
description of the dynamical state of the RRE. Such a timely dynamic information can be obtained 
by using peripheral Application Specific Processors (ASP) interacting with a general-purpose Central 
Processing Unit (CPU), i.e. the expert system, in the spirit of sixth generation computers [2]. The 
ASP would include hardwired, parallel distributed architectures, e.g. neural networks, to extract 
global dynamic state information from a myriad of sensor data (in the time and space domain). 

This cybernetics perspective of the ICS is suggestively illustrated by analyzing the visual infor- 
mation processing pathway: eye-ball/retina/brain. The sensor function is performed by a spatial 
filter (eye-ball), which focuses the individual photons onto the retina (ASP), where global informa- 
tion such as contrast enhancement and motion detection is extracted by an intricate neural network 
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of cones, rods, bipolar and horizontal cells. When reaching the brain, (general purpose CPU), the 
visual information has been compressed in a ratio of 150/1, and therefore can be more efficiently ana- 
lyzed in relation to patterns or information previously stored in the brain. Likewise, the engineering 
model stored in the expert system can be best utilized on the premise of an efficient compression of 
the sensor data. 

From this perspective, we propose an integrated architecture which will regularly measure the 
load values as a function of time, identify the presence of load cycles, and extract in real-time their 
contribution to the cumulative fatigue damage. Component life prognosis and failure detection will 
be achieved through such an accurate monitoring of the load history. 


2 Architecture of the Real Time Fatigue Life Estimator. 

As indicated in Fig. 4, the real time fatigue life estimator that is proposed consists of two functional 
blocks: a preprocessor and a neural network. The role of the preprocessor is to perform simple 
front-end signal processing in order to improve the computational efficiency of the neural network. 


2.1 Preprocessor. 

First, the preprocessor extracts the peak-to-peak load transitions by comparing the loads sampled 
at three successive time intervals. Once an extremal load value has been identified, it is stored in 
a shift-register type of buffer which is continuously analyzed to detect the occurence of load cycles. 
If a load cycle is identified, the content of the shift-register is input to a neural net architecture. 
If not, the shift-register is ready to receive an additional extremal load value. In order to define a 
load cycle and understand its impact on the fatigue life, it is best to review briefly the underlying 
concepts of the uniaxial local strain approach [3] [4] [5]. 

2.1.1 Uniaxial Local Strain Approach. 

To demonstrate the concepts of the local strain approach, it is suggestive to interpret the mechanical 
effects of load changes in terms of dynamical trajectories in the stress /strain plane. It assumes a 
uniaxial behavior of the component, i.e. a stress cr and a strain e, in the direction of the applied 
load P. Changes in stress are related to changes in strain according to the standard elastic/plastic 
decomposition of the total strain range Ac 


A ^<7 

A< = ~E + 


( 1 ) 


where the elastic strain range, , is the ratio of the total stress range, Act, over the elastic modulus, 
E i the plastic strain range. This non-linearity of the stress/strain relationship is responsible 

for the occurence of hysteresis loops in the {a, e} plane which characterize the degree of fatigue 
damage of the component caused by the corresponding load variations. Examples of trajectories in 
the {cr, e} plane are given in Fig.l for an irregular load sequence applied to the fatigue specimen of 
Ref[3]. Properties of RQC100 steel are used in this paper to illustrate the correlations between load 
sequences and fatigue damage. A peak value of the load corresponds to a reversal in the loading 
direction, e.g. compression to tension, or vice-versa. As a result of such reversals, small cracks are 
initiated and enlarged, leading ultimately to mechanical rupture. The fatigue damage of a load cycle 
is therefore characterized by the number of times, N cyc i ei that such a cycle can be repeatedly applied 
to a component before a crack develops in the material. For this reason, we will refer to A ; as 
cycle number. 
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According to the Linear Damage Rule (LDR) [6]-[7], the fatigue damage p cyclt of a hysteresis 
loop (such as one of those in Fig. lb) is given by the ratio 

1 


1 ^ Pcyclc — 


IV, 


cycle 


> o 


( 2 ) 


An estimation of the fatigue life is obtained by cumulating the contributions p cycle of the indivi ua 
cycles which occured during the load history. In the LDR approach, the fatigue damage of a load 
cycle is independent of the fatigue damage accumulated by the component. This cycle-coun ing 
method is useful for the treatment of random load sequences, where the load history consists of a 
succession of cycles randomly distributed in the stress/strain plane. Yet, its predictions are non- 
conservative when the interaction of loadings at different strain/stress leve s has to be aken into 
account. Such is the case of high-low load sequences, where the fatigue damage of a load cycle 
not only depends on its cycle number , N cycle , but also on the * ** 

component. This necessitates the use of Damage Curve Approaches (DC A) [8]- [9] to reflect the 
effects of such non-linear interactions on component fatigue life. In spite of these non-linearities, 
component fatigue life can be estimated from the sole knowledge of the mapping between load cycle 
and cycle number, whether the Fatigue Damage Accumulator is implemented using LDR or DCA 

^IUs the purpose of this paper to demonstrate the possibility of learning, based on a local strain 
approach, the mapping of random load cycles with their fatigue damage p cyc ie ^cycl t 

The way the preprocessor detects the existence of a load cycle from a set of four consecutive 
extremal values is indicated in Fig.2, together with a diagram of the corresponding hysteresis loop 
in the {o-,c} plane. In Fig. 2a, the compression sequence AB is interrupted in B by the hysteresis 
loop BCB (tension sequence BC, followed by compression sequence CB), before resuming to BD. 
Fig. 2b is obtained from Fig. 2a by exchanging compression and tension. The existence conditions 
for hysteresis loops are P A > Pc and P B > Pd (Fig.2a), or P A < Pc and P B < Pd ( Fl S- 2b )' 
Compressive and tensile load cycles are shown in Fig. 3a for a load sequence resulting rom e 
superposition of low-frequency and high-frequency vibrations. In contrast, several load cycles may 
be initiated without being completed, as is the case for the damped vibrations shown in Fig. 3b. 

Whenever load reversals associated with a load cycle have been identified, the preprocessor passes 
the result of the neural computation to the Fatigue Damage Accumulator as shown in Fig. 4; and 

then removes the load cycle from the content of the shift-register. _ 

Although the local strain approach assumes the existence of a uniaxial state of stress, its un- 
derlying concepts of hysteresis looping and memory effects are sufficiently general to be valid under 
more complex conditions [3], e.g. biaxial/triaxial states of strain/stress. 


2.2 Neural Network Architecture. 

Once a load cycle has been identified by the preprocessor, the contents of the shift-register are 
input to a feedforward neural net. As shown in Fig.4, each input unit of the neural net is in 
one-to-one correspondance with each buffer of the shift-register of the preprocessor. This means 
that a peak value of the load sequence is fanned out to each neuron of the first hidden layer. 
The activation function of a neuron is chosen to be the exponential type of I/O response function 

output = 1 — where the input is the sum of the weighted outputs of all the connected 

neurons and IT threshold term. The neural net output contributes to the fatigue damage only if the 

presence of a load-cycle is identified by the preprocessor. 

An important issue which affects the design of the neural architecture is the accuracy with 
which the fatigue damage needs to be estimated. To estimate the cumulative fatigue damage o a 
component, it is necessary to estimate accurately the fatigue damage due to each individual cycle 
since the load history consists of a multitude of load reversals of small amplitude. We propose 
a neural architecture' that can be expanded laterally in order to adaptively increase the precision 
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of the neural computation. The architecture is an arrangement of neural net building blocks that 
is functionally reminiscent of a successive approximation register, operating in parallel and it is 
described as follows. 

If ' Pmin is an acceptable lower bound for the fatigue damage of an individual hysteresis loop a 
eed-forward net Fi is trained to map individual load cycles over the interval [0, -log(p mtn )l. If F 
has learned the mapping with a sufficient accuracy, it can be used to predict the ’fatigue"life. If no/ 
Fi is used as a classifier to separate the cycle contributions that are larger than -0.51og(p mtn ) from 
those that are smaller. Two feed-forward nets, F 2 and F 3 , are then trained separately to map the 
load cycles over the intervals [ 0 , + 6} and ^ - log{pmtn)] respectivelVi where 

6 is the average resolution of F,. After training, F u F 2 , and F 3 are operated in parallel, with the 
output no e of F\ inhibiting the output node of Fj or F 3 which does not correspond to the class 
defined by F, As shown in Figure 5, the remaining active node of F 2 , or F 3 , provides an estimation 
of the fatigue damage caused by the load-cycle which has been identified by the preprocessor If this 
estimation is accurate, the parallel arrangement of F u F 2 , F 3 can be used for prognosis/diagnosis of 
the mechanical fatigue. If not, F 2 and F 3 are used as classifiers, and the architecture is expanded 
laterally with four additional nets whose outputs are activated/inhibited by the output nodes of F 2 
^3' This procedure can be reiterated until desired accuracy is reached. 

To illustrate the practilality of this approach and demonstrate the possibility of extracting the 
highly non-linear features which govern the material behavior in complex situations, an architecture 
of two-hidden layer neural building blocks has been trained in supervised learning mode [lOl-flll 
*° pre< J lct the fati S ue damage of a specimen of RQC 100 steel in the uniaxial local strain approach. 
When learning from experimental data, the choice between a back-error-propagation network (BEP) 
or a counter-propagation network will depend upon the degree of representativeness of the training 
data set. When learning from an analytical model, e.g. the local strain approach, it is possible to 
generate a training data set that is highly representative of the load/damage relationship. In that 
case, a good internal representation of the mapping can be developed by a BEP network [12]-[13] 
known to be more efficient as the degree of representativeness of the training data set increases. 


3 Neural Network Training. 

The training data set was generated on the basis of the uniaxial local strain approach for cyclic 
loading within fixed limits [3] [5]. The initial value of strain, e 0 , is calculated according to the load 
strain calibration of the loading path 


e 0 = sign(P 0 ) x (J^! + (^) ' ) ( 3 ) 

where P 0 is the final value of the loading path, and the fitting constants are Cj = 1410 ksi (kilopounds 
per square inch), C 2 = 65.5 ksi, and d = 0.31 for the geometry of the uniaxial RQC 100 specimen 
given m Ref.[3J. The initial value of the stress a 0 corresponding to F 0 is calculated following the 
non-linear stress/strain relationship given in Eq.(l^) 




X/s 


sign(a 0 ) = si 0 n(€ O ) ■ 


(4) 


where the elastic modulus is E = 29500ksi, A = 167ksi and s = 0.10 for RQC 100 steel. The initial 
loading creates a cnstallographic slip in the material, which modifies its mechanical response under 
su sequent loadings. The irreversible effect of this cristallographic slip on component behavior is 
accounted for by expanding the calibration curve and the stress/strain curve, Eqs.( 3 )-( 4 ) with a 
factor of 2 and by shifting their origins. 
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After loading, the ” image” of a load sequence in the {<7, e} plane is therefore obtained as the 
solution of the set of equations 


|g - e r | _ 1 -P - ?t\ A P - Pr\ , 1/d _ k - grl , \f_ ~ fV 

2 2Ci 1 2C 2 2E [ 2A 


sign(a — a r ) = sign(e — e r ) = sign(P — P T ) 


(5) 


where (P r ,e r ,cr r ) corresponds to the origin of the stress-strain hysteresis loop curve of the current 
’’point” (P, e,<r). For a hysteresis loop #12 starting at (Pi, € 1 ,^ 1 ) and ending at (P2, ^2^2)} the 
number JV 12 of loops # 12 that would lead to the initiation of a crack is obtained from the fitted 
equation [3] 


l f 2 ~ f ll 
2 


(2N l2 ) b + ei f (2N 12 ) c 


( 6 ) 


with cr m = <T : * <7 2 if it is a tensile mean-stress, he. > 0, and cr m = 0 if it is a compressive mean- 
stress. For RQC100 steel, the fatigue strength coefficient is atj = 168ksi, and the fatigue strength 
exponent is b — —0.075; the fatigue ductility coefficient is elf = L06, and the fatigue ductility 
exponent is c = —0.75. 

It is well-known that the performance of the neural computation can be significantly enhanced 
by embedding explicitly any a-priori knowledge or constraints of the system in the preprocessing or 
the neural architecture itself. The proposed neural estimation of the fatigue life is actually subject 
to the detection of load cycles by the preprocessor. Similarly, any a-priori knowledge about the gross 
dynamic properties of the system is likely to improve the performance of the neural computation. For 
example, load sequences that are periodic or quasi-periodic, such as in Fig. 3a, lead to a maximum 
of two or three cycles that remain simultaneously open during the entire load history. This implies 
that the mechanical behavior of components subjected to load patterns such as the one of Fig. 3a 
can be predicted by reducing the number of input units to a maximum of 8-10. In addition, the 
regularity of such load patterns facilitates the feature extraction, and improves the ability of the net 
to predict component behavior for such load patterns. In contrast, for the most irregular sequences 
where the loads are randomly distributed, load patterns such as in Fig. 3b can occur, so that more 
buffers are needed in the shift-register of the preprocessor. In the latter case, training/learning is 
more difficult because of the great irregularities of the load variations. 

The behavior of RQC100 steel under random load sequences has been simulated for cyclic loading 
within fixed limits [3] [4] [5] in the dynamic range D p = [— Pmax , *F Pmax] • The initial load Pi is 
arbitrarily chosen over D V) and defines the domain of variations of the subsequent loads P(i), 
|P(t)| < |P»|. (The neural estimation of cumulative damage resulting from cyclic loading with 
variable boundaries, for which there exists t such that |P(t)| > |P<|, is the object of an on-going 
analysis.) The random load sequences used to generate the training data set are obtained as follows. 
A uniform set {Pi} of initial loading values is generated over Dp. For each P t , a uniform set {P a v} of 
average values of the load sequence is generated with the condition that \P{\ > P av > — |P*|. Given 
an average value P av , a uniform set {P a } of maximum amplitudes of the load variations is generated 
over [0, max(\Pi -b P av |, |P* — P at ,|)j. For each triplet (P t , P atn P a ), a load sequence is generated as a 
set of K points randomly distributed over [P at; — P a , P ar + P a ] , from where only individual load cycles 
such that 1 > pcyde > pram are retained for training. The interval to be mapped, [0, —log(p m i n )], 
is divided in 10 equal subintervals {Ik.k = 1, 10} which are equally populated by data generated 
from random load sequences. A training iteration consists of choosing randomly a value —log(p cyc i e ) 
in one of the randomly picked subintervals /*• The load values of the corresponding sequence are 
scaled as P cy cic = p cyctc e[— 1. 4-1], and the fatigue damage is scaled over the active range [0.1, 0.9] 
of the neuron output of the last layer: 


-target 
P cycle 


x 0.8 + 0.1 

log{pmin) 


(7) 


5 



To favor homogeneous learning of the mapping load/damage over [0, — log(p rnin )] 1 the weights and 
thresholds were updated by batching the steepest-descent contributions of 10 consecutive I/O pairs 
(Pcycle,p; v T)- 

In this simulation of the behavior of a uniaxial RQC100 steel specimen [3], the dynamic range 
of load variations has been bounded by P ma x = 90ksi, i.e. , approximately at half the value of the 
fatigue strength coefficient <j\ j — 167ksi. With this value of jP max , it was found that load sequences 
corresponding to K = 200 random values were long enough to indicate that a probable maximum 
number of 7 initiated cycles can remain simultaneously open during a random load history. On 
the basis of this statistical result, 15 buffers were allocated to the shift-register of the preprocessor. 
Cycles contributing less than p min = 10“ 10 were neglected, and 10,000 elements were used for 
training. 

4 Neural Network Performance. 

Training was done on a PC-based accelerator coprocessor board (ANZA-f) and the elementary 
neural building block was a two-layer feedforward net with 15 input units, 100 neurons in the first 
hidden layer, 50 neurons in the second hidden layer, and a single neuron in the output layer. The 
performance of the neural computation was evaluated by averaging the absolute value of the error 
tcycie = | P'cyciT — I over a set ° f 10,000 I/O pairs (Pcyciei Pc^cie 1 ) that were not used to train 

the net, and with the scaled fatigue damage uniformly distributed over [0.1, 0.9]. 

The average absolute error of this net Fi is e —< e cyc ie > — 0.021 with a standard deviation 
6e = 0.029, (versus the average ei tarnxng = 0.020 and the standard deviation 6ei ea Tnxng = 0.026 
calculated for the training data set); and the distribution of the estimated values is symmetric 
with respect to the exact values. This represents accurate learning (within 2%) of the uniaxial 
strain approach. To interpret the effect of this error on fatigue damage estimation, recall that a 
logarithmic relationship exists between the network output and the fatigue damage estimate. Thus, 
the fatigue damage predicted by the network has the same probability of being larger or smaller 
than that predicted by the teacher. If it is larger, it will be in the average r ~ 1.8 times the 
value predicted by the teacher; and if it is smaller, it will be in the average ~ times the 

value predicted by the teacher. Said another way, 90% of the estimated values will be within a 
factor 4 and 1/4 (standard deviation upper-bound) of the exact data. In order to demonstrate the 
ability of the architecture to improve the accuracy of the neural computation, two additional neural 
building blocks, F 2 and F 3 , have been trained to map load-cycle/fatigue-damage for [0,5.5] 

and 10] respectively. As shown in Fig. 5, net F\ is used as a classifier for the output nodes 

of F 2 and F 3 . With this two-stage neural architecture, the value of the average ratio r defined above 
is lowered to r = 1.5, and 90% of the estimated values are within a factor 2.5 and 1/2.5 (standard 
deviation upper-bound) of the exact data. It is also noted that 75% of the estimated values are 
already within a factor 1.7 and 1/1.7 of the exact data. 

It is finally emphasized that no attempt has been made to look for the minimum number of 
neurons and synaptic connections that provide the same degree of accuracy as the two-hidden layer 
building block chosen in this simulation. 

5 Conclusion. 

A neural net architecture has been proposed to extract the correlations between random load se- 
quences and their contributions to cumulative damage in order to predict component fatigue damage. 
The efficiency of the neural computation to map load sequences with fatigue damage has been demon- 
strated on a uniaxial RQC100 notched member for the most irregular load histories of cyclic loading 
with fixed boundaries. The real-time capability of this parallel processing architecture, together with 
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the simplicity of its practical implementation, makes it computationally very attractive to enhance 
P" formance of com P lex ICS. Because of the complexity of the component couplings of the RRE 
the ICS has to be able to anticipate in real-time the effect of global changes of conditions on indi- 
vidual components, e.g. the effect of a gas flow variation on the fatigue life of a turbopump blade 
In that context it would be natural to extend the concepts of neural computation developed in 
MAI hci thC fatlgUC hfe P redlctlon of muUiazial components subjected to temperature variations 
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(a) Loads distributed betweeen *10 and 10 kips. Initial 
loading value of 10 kips. 



STRAIN, t 

(b) Note the closed hysteresis loops 5-8-5’ and 6-7-6\ and the 
gradual formation of the outer cycles. 

Figure 1 . - Open load sequence and its image in the {o,e} plane. 
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(a) Compressive. 




(b) Tensile. 

Figure 2. - Load cycles and their image in the (o, e] plane. 




(b) Irreducible load squence. 

Figure 3. - Extreme cases of load sequences. 
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